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We investigate the full 5D dynamics of general braneworld models. Without making 
any further assumptions we show that cyclic behavior can arise naturally in a fraction of 
. . . physically accepted solutions. The model does not require brane collisions, which in the 

■ stationary case remain fixed, and cyclicity takes place on the branes. We indicate that the 
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cosmological constants play the central role for the realization of cyclic solutions and we show 



^ i' that its extremely small value on the observable universe makes the period of the cycles and 

, the maximum scale factor astronomically large. 

in 
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I. INTRODUCTION 

p) The last decade proves to be really exciting for cosmology. Observational data indicated, among 



cd: n 

, other very interesting results, that the expansion of the universe is accelerated At the same 

time the braneworld scenario appeared in the literature M, yl . Though the exciting idea that we 
CN ■ 

i live in a fundamentally higher-dimensional spacetime which is greatly curved by vacuum energy 

' was older {4], the new class of "warped" geometries offered a simple way of localizing the low energy 



gravitons on the brane. 

n this novel background the old idea of a cyclic Universe was reheated. Started as ekpyrotic 



T— I in tnis novel bacKground tne old idea 01 a cyclic Umvers 

^ : fl 1 enriched to ekpyrotic/cyclic Q, y, H Q, Q, Q 



and recently to new ekpyrotic 



14. jla. [1^ . the new paradigm tries to be established as an alternative to standard cosmology. 
According to its basic contents, our universe experiences an infinite or extremely large number of 
cycles, each one consisting of a hot bang phase, a phase of accelerated expansion, a phase of slow- 
ekpyrotic contraction and a bounce-bang that triggers the next cycle. Starting with a simplified 
notional framework (infinite and not "created" time) cyclic cosmology have many advantages. It 
successfully faces the homogeneity, isotropy, topological and flatness problems, it handles the issue 
of initial conditions, it incorporates the dark energy and transforms it to an important factor, and 
it provides the mechanism of the generation of cosmic perturbations and of structure formation. 
However, there are some key issues that do not have a consistent and efficient approach so far, 
despite the great progress. These are the settlement of the singularity, although temperature and 
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density remain finite, the entropy evolution, and the fate of the perturbations through the bounce. 
Through this research, cychc scenarios have become more comphcated, by the insertion of more 



complex potentials, of more branes of the mechanism of ghost condensation [13I. \y\. of 



more 



scalar fields 



14( 1 and of procedures which cancel the tachyonic instabilities |l5l ]. 
Most of the works on cyclic cosmology involve, initially or at some stage, the transition to 
effective AD equations. However, as it was mentioned in 13], such a procedure does not lead 
to reliable results since one cannot return to the 51? description self-consistently. Furthermore, the 
old 4L'-singularity problem (of both Big Bang and traditional cyclic universes), has been replaced 
by a new one (equally annoying) concerning the singularity of extra dimension(s). This later case 
is accompanied by the brane collision phenomenon, which seems to be a basic constituent of the 
ekpyrotic scenario. 

In this work we desire to investigate the full 5D dynamics of general braneworld models and 
examine if a cyclic behavior is possible. This is an essential procedure in order to consistently 
confront the arguments of the authors of [isl, which claim that cyclic behavior cannot arise from 
a complete 5D description, and our study must not include any additional assumptions or fine 
tunings in order to remain general and therefore convincing. Secondly, we are interested to explore 
if a cyclic behavior of 5D dynamics is necessarily related to brane collisions. This work is organized 
as follows: In sectionlUwe present the 5D braneworld model and we derive the equations of motion. 
In section Hm we provide analytical solutions for two simplified stationary solution subclasses, while 
in IIVI we investigate numerically the full stationary dynamics. Finally, in section |V] we discuss the 
physical implications of our analysis and we summarize the obtained results. 



II. THE MODEL 



We consider quite general braneworld models, characterized by the action 



2l|: 



d'^xdy^ 



gR+ I d xdyy/^ 



1(50)2 -y(<A) 



[ d'xv^{[K]+u.m, (1) 

1=1,2 "^^^ 



where = -jj^ is a 5D gravitational constant, and all quantities are measured in units of M5. The 
first term describes gravity in the five dimensional bulk space. The second term corresponds to 
a minimally coupled bulk scalar field with the potential V{(j)). The last term corresponds to two 
3 + 1 dimensional branes, which constitutes the boundary of the 5D space. We allow for a potential 
term U{(j)) for the scalar field at each of the two branes, and we denote by 7 the induced metric on 
them and by K their extrinsic curvature. Here and in the following the square brackets denote the 
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jump of any quantity across a brane {[Q] = Q{y+) — Q{y-)). The reason we use a second brane 
is to eliminate possible "naked" singularities, therefore by assuming /Z2 symmetry across each 
brane we restrict our interest only in the interbrane space. 

As usual the two branes are taken parallel, y denotes the coordinate transverse to them and we 
assume isometry along three dimensional x slices including the branes. For the metric we choose 



the conformal gauge 



2i: 



ds^ = e2^(*'^) {-de + dy^) + e2^(*'^)(ix2. (2) 



2C 



This metric choice, along with the residual gauge freedom (t, y) — > {t',y') which preserves the 
2D conformal form, allows us to "fix" the positions of the branes. Without loss of generality we 
can locate them at y = 0, 1, having in mind that their phy sical distance is encoded in the metric 
component B{t,y), and at a specific time it is given by 



D{t)^ / dy^^= / dye^^''y\ (3) 
Jo Jo 

quantity that is invariant under the residual gauge freedom in our coordinates. The reason we 
prefer the metric ([2]), instead of the usual form in the literature, is that in the later case the brane 
positions are in general time-dependent and the various boundary conditions are significantly more 
complicated. Thus, our coordinates are preferable for numerical calculations, despite the loss of 
simplicity in the definitions of some quantities. Eventually, the physical interpretation of the results 
is independent of the coordinate choice. 

The non-trivial five-dimensional Einstein equations consist of three dynamical: 

A-A" + 3i2 _ 3A'2 = -e^^y 

3 

B-B"- 3i2 + 3A'2 = + ^ - le^By 



2 2 3 

^ - + 2,A4) - ^A'cj)' + e^^V^ = 0, (4) 

and two constraint equations: 

- A' A + B'A + A'B - i' = -(/.(/)' 

3 

2^'2 _ ^/^/ + _ ^2 _ ^ _^ _ ^ _ \ 2By .5) 

6 6 3' ^ ' 

where primes and dots denote derivatives with respect to y and t respectively. It is easy to show 
that the constraints are preserved by the dynamical equations. 
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Additionally, from the boundary terms in the action for the branes we obtain the following 
junctions (Israel) conditions: 

m = ^1 Ue^ 

[</.'] = ±e^C/,<^, (6) 

where the upper and lower signs refer to the branes at y = 0, 1, respectively. Since we have imposed 
Z2 symmetry across the two branes, for any function Q we get: 

[Q']o = 2Q'(0+) , = -2Q'(1-). (7) 

For the bulk potential V{(f)) we assume a general form consistent with the stabilization mecha- 
nism {221: 

= ^mV + A, (8) 

where A is the 5D bulk cosmological constant. For the brane potentials Ui{(j)) we use the quite 
general quadratic form 



21|: 



Ui{^) = ^Mi{<l), - a,)^ + Xi, (9) 

where Aj stand for the brane tensions and (pi for the value of (p on the i-th brane. The "masses" m 
and Mi can be varying. 

Finally, the induced 41) metrics of the two ("fixed" -position) branes in the conformal gauge are 
simply given by 

ds'^ = -dr^ + a\T)dx.'^, (10) 

with dTi = e^^dt and Oj = e^^ the proper times and scale factors of the two branes. Thus, for the 
Hubble parameter on the branes we acquire 

1 da 



a dr 



Ai, (11) 



which is invariant under residual gauge transformations. As usual, we identify the brane at y = 
as the visible brane corresponding to our Universe. Note that relations (jlOp and (jlip can be 
generalized to hold in every 3-1-1 slice transverse to the fifth dimension. 

Finally, let as comment on the behavior of gravity on the physical brane in our model. As 

n 

was shown in [23], in such a two-brane model we re-obtain the correct Newton's law in the brane- 
universe. This becomes more transparent if we include the aforementioned brane stabilization 
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mechanism 22], where the two-brane model leads naturally to the recovery of the Einstein gravity 
on the physical brane [2^, [2^ • More generally, even in a generic two-brane model with arbitrary 
bulk and brane potentials, using the renormalization group flow of the 4D Newton's constant in IR 
and UV, it can be shown that Newton's law is recovered, plus one extremely small brane correction 



25| 



III. STATIONARY CASE. ANALYTICAL SOLUTIONS 

The model we have described is general enough and includes the full spacetime evolution of 
the 5D braneworld. Our aim is to look for physically accepted solutions that correspond to cyclic 
behavior. The main difficulty of solving the equation system ([4])-(l9]) is that one has to satisfy the 
constraints ([5]) and boundary conditions ^ at t = to, ensuring that no divergencies are present in 
the interbrane space. Then, one has to assure that time evolution will not be unstable, or give rise 
to naked singularities in the bulk. 

Let us first investigate a subclass of solutions, the so-called stationary ones. In this case we 
assume that 

A{t,y)^B{y)+Ht. (12) 

This case is characterized by a fixed bulk geometry and maximally symmetric (de Sitter, Minkowski 
or oscillatory) branes: 

ds^ = e2^(f) {dy^ - dt" + e2^*dx2) . (13) 

Note that there is a name confusion in the literature, since spacetimes (I12p are called stationary 
i n 26( 1 and when H = they are called static (a convention that we follow in this paper), while in 
20] they are called static even for H ^ 0. 

Under (fT2l) the system of equations (jl]) and ([5]) transits to the time-independent form: 

B"iy) = -^-e'^iy) ^"^y^ 



6 4 
4>"{y) + W'{y)4>'{y)-e'^(y^^ = 

k'2 



H' = B'\y) + \e^''^y)V[<t>)-^, (14) 

where (p is also assumed to be time independent. Furthermore, boundary conditions ([B]) will be 
satisfied at all times, provided that they do so at t = tg- Finally, according to definition (jll|) and 
under the ansatz (|12p . the Hubble parameter on the physical brane is simply Hq = e~'^^^^^ H'^ . 
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Equations (jl4p and boundary conditions ^ cannot be solved analytically for the general bulk 
and brane potentials of ([8]) and ([9]). However, this is possible in two simplified cases: A) Assuming 
V{4>) = with full Ui{(t>), and B) assuming V{(j)) = A, Uo{(j)) = Aq and Ui{(j)) = Ai. In the following 
subsections we examine these two cases successively. 



A. V{(f>) = and U^{(p) = \M, 



-(7iY+ A, 



If we set V{cl)) = we can acquire analytical solutions depending on the sign of H^. For > 
and setting H = \ V H'^\ we acquire: 



B{y) = B{0) + ^ log sinhSHy + coshSi^y 



and 



0(y) = 0(O)-^log 



{e^^y -u){l + u) 



(15) 



(16) 



where we set u = \/ i Biud |i?'(0)| > H always in this case if we want (/''(O) to be real. 



For < and setting 6 = \\/—H'^\ we get: 



B{y) = B{0) + ^ log (^-^ sin SOy + cos 39y 



(17) 



and 



<Piy) = m-^^og 



tan ( ^6y + i arctan W 



B'(0) 



(18) 



Let us make some comments here. Firstly, in the expressions above we have eliminated the 
three integration constants in terms of B(0), B'{0) and 0(0), since it is more convenient, but one 
could equally use the values of three of B, B', (p and (p' at any y. Moreover, these solutions have 
to satisfy the boundary conditions ([6]). Eventually, all quantities, including H^, are given in terms 
of the six parameters of the model (Mq, Aq, ctq. Mi, Ai, cti, since we have set m = and A = 0). We 
will return to this subject in the next section. Secondly, note that there is analytical continuation 
between the solutions for positive and negative H^, i.e. the two solutions coincide for H — > i9. 
Thirdly, it is easy to see that in the limit — > 0, B(y) and (j){y) become smoothly a negative 
logarithm of a linear in y function. Lastly, note that due to the non-linear form of the differential 
equations there exist more solution branches which can be obtained straightforwardly from ()15|) - 
(llSp by sign changing. However, since all of them imply the same metric for the visible brane, 
which is given below and which is the central subject of this work, we do not write them explicitly. 
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In order for the solutions to be physically accepted, we have to assure that no naked singularities 
are present in the interbrane space, i.e. between y = and y = 1. The investigation of the 



> case was presented in 27|] where we showed that values are restricted to very small 
values and only one 5D sub-surface of the 6D parameter space allows for arbitrary large 
values. Here we repeat the cogitations in the < case. Solutions (fTTl) and (fTH]) have poles at 
sin 39yp + cos 36yp = 0, i.e. at 



with n G Z, where yo = arctan y~-gTjQ^j with —it < yo < 0. In this notation ypo stands for the 
largest negative pole. Demanding none of these ypn lying in [0,1] we require that the smallest 
positive pole, i.e. y^i, to be greater than 1. This leads to yo > 30 — tt and since — vr < yo < we 
acquire: 



361 - vr < arctan y-^rjjyjj < 0- (20) 

Thus, this relation provides a narrow and absolute window for 6, in a sense that there are no areas 
at all in the 6D parameter space that give 9 larger than vr/S. 

Finally, we desire to explore the forms of the bulk and brane geometries in these cases. For 
> we get de Sitter, anti-de Sitter or Minkowski branes (3+1 slices in general), that have been 
studied extensively in the literature. The bulk structure is fixed and this can be also confirmed 
by calculating the brane physical distance according to ([3]), which is obviously time independent 
(given analytically as an expression in terms of hypergeometric and hyperbolic functions of the 
parameters), i.e. the branes do not move. 

For < 0, which is the case of interest in this work, we acquire also a fixed bulk with constant 
brane distance. For the branes we get an interesting oscillatory behavior. In particular, the 5D 
metric: ds'^ = e^^^^^ i^—dt^ + dy'^ + e^*^*(ix^), implies a AD metric of the form (jlOp on the visible 
brane, with the scale factor given as: 



a^ir) = c sin^ 



(e-«(°)0r) , (21) 



where r is the visible brane proper time defined above as dr = e^^'^^dt (for clarity in the following 
we omit the index i which distinguishes the proper times of the two branes, using r for the visible 
brane and ri for the hidden one). The constant c depends on the exact relation between r and bulk 
time t (it is equal to e^^^^^ for r = e^^^H). Moreover, one should add a constant phase which would 
be determined by the a^(T) value at a specific r. Without loss of generality we have performed 
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a shift in r such that a^(0) = (note that r = is not an initial time but just a random one). 
Definitely, one could achieve a relation similar to ()21|) in terms of bulk time t, but we desire to 
present the results from the visible brane observer's point of view. 



B. V{cl>) = A, Uoicb) = Ao and Ui{^) = Ai 



In the case where the branes have just tensions, the bulk potential is comprised of a non-zero 
cosmological constant and there is no scalar field, analytical solutions, depending on the sign of 
H^, can be derived. 

For > 0, with H = \ V we acquire: 



B{y) = -log i + 



2H^ 



12ci 



(22) 



where ci = 



'B(O) 



23]. Expression ([22]) . 



similar to the results of 
apart from y ±00, has a naked singularity at yp = 277- log (—-^^) which must lie outside 
[0, 1]. It is obvious that when A < it is singularity- free and is the well studied AdS bulk case 
of the literature. Furthermore, solution (j22p . in the specific case where A < 0, Aq = \/— 6A and 
Ai = —\/—6A, corresponds to the two-brane Randall-Sundrum model j^, where acquires a zero 
value and B{y) becomes a negative logarithm of a term linear in y. 
For < 0, and setting 9 = \ V^IP\ we get: 



Biy) = log 



6 sin 9y + arcsin 



(23) 



Note that in this case A must be negative, as it is easily implied from the last of equations ()14p if 
we demand B{y) to be real, i.e. the bulk is always AdS. Expression (|23p possesses singularities at 



(24) 



with n G Z, where yo 



arcsin 




with — TT < yo < 0. As previously, ypo stands for the 



largest negative pole. Forcing the smallest positive pole, i.e. ypi, to be greater than 1 we acquire 
yo > 9 — TT, and since —n < yo < we finally obtain: 



TT < — arcsm 



< 0. 



(25) 
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Thus, in this case there are no areas at all in the 3-D parameter space that lead to solutions with 
9 larger than vr. 

Let as describe the spacetime geometry that corresponds to these solutions (we repeat that 
there are more solution branches arising form (j22|) and (j23p by sign changing, but they correspond 
to the same brane metric). In the > case we acquire two stabilized branes, with their 
time-independent physical distance given analytically through ^ in terms of inverse hyperbolic 
trigonometric functions of the parameters. The AD induced geometry on the two branes can be 
de Sitter, anti-de Sitter or Minkowski and has been investigated widely in the literature. For 

< we also obtain a fixed bulk with the time-independent physical distance of the two branes 
given easily analytically. However, the geometry on the branes acquires the interesting oscillatory 
behavior described in the previous subsection, with a scale factor of the form: 

a\T) = c sin^ (e-^W^r) , (26) 

identically to relation (j2ip . 

Relations (j2ip and (j26p correspond to cyclic Universes. Indeed, they imply that the scale factor 
pulsates between the following extremum values: 

amax = Vc, (27) 

where as we mentioned c depends on the specific relation between t and r (c = e^^^^^ for r = e^^^H). 
The constant period of the cycles is given as: 

r = ^. (28) 

In (|27p . (j28p all quantities are measured in units of M5. The aforementioned behavior corresponds 
to a sequence of expanding and contracting phases on the branes, while the physical brane positions 
stay constant. The fifth dimension remains unaffected, that is there are no brane collisions. The 
bulk dynamics determines also the quantitative characteristics of the cycles. Indeed, the value of 
B{y) on the physical brane specifies the maximum value of the scale factor and can be arbitrary. 
Moreover, it designates the cycles period. Fortunately, the fact that 9 is bounded from above 
(according to (j20p and (j25p l restricts completely its effect on decreasing the period. We will return 
to these subjects in the discussion section below, where we show that only astronomically large 
maximum scale factors and periods are possible. Finally, note that in the aforementioned analysis 
cyclicity arises naturally from the 5D dynamics, and the oscillatory behavior of the scale factor is 
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explicit and not a result of Hubble constant sign change in special cases 



flln 



13|, 



291 ]. Furthermore, 



301 ] or the assumption 



it emerges from general braneworld models, without Gauss-Bonnet terms 
of charged AdS bulk black holes which charge (along with fine tuning) is responsible for restricted 
cyclic behavior [lo| . 

We close this section by referring to the stability of the aforementioned solutions in both metric 
and scalar field components. As it was shown in l3ll |. a two-brane model with the bulk and 



brane potentials considered above, has stable solutions, free of tachyonic modes. In particular, the 
presence of the stabilization mechanism |2^, which forbids branes to move, makes all physically 
accepted solutions (i.e. those which are singularity-free and satisfy boundary conditions) to be 



24 ]. This feature has been also confirmed by numerical investigation 



stable under perturbations 
□ 

in [20(] and has been verified by us, too. Thus, the stability of the solutions provides the necessary 



physical hypostasis to our model. 



IV. STATIONARY CASE. NUMERICAL RESULTS 



In the previous section we derived analytical solutions for the full 5D equations, in the stationary 
ansatz with two simplified potential cases. We expressed the solutions in terms of the values of 
-B(O), B'(0) and (/){0) and we stated that these must satisfy the boundary conditions Q. The 
usual, in the literature, method to achieve this is to randomly choose and the solution values 
at y = and then fine-tune the model parameters. However, this procedure definitely does not 
reveal the properties and the rich structure of the solutions. On the contrary, as we showed in 
it restricts the investigation in a small subclass of solutions. The natural way to encounter the 
problem is to choose randomly and uniformly the potential parameters at first, and then seek for 
physically accepted solutions, i.e. divergencies-free expressions which satisfy equations (jl4p and 
boundary conditions ([6]). This is a hard task in general. The method we use is the following: We 
first choose randomly the values of the model parameters, uniformly distributed in a hyper-cube 
(6D in the case of subsection IIII Al and 3D in that of IIIIBp . The obtained results do not depend 
on the hyper-cube's size, but on its effectual covering (number of parameter multiplets used in 
the calculation). We use the globally convergent Schmelcher-Diakonos algorithm [32] to solve the 
transcendental equation system with accuracy 10^^^. We find hexads (or triads) of parameters 
corresponding to acceptable solutions, and calculate Hq through Hq = e^'^^^^^ H'^ . 

Numerical investigation reveals that only a small fraction of parameter choices (we use 10^ 
parameter multipltes) allows for solutions to exist (~ 15% in case IIII X\ and ~ 5% in case IIIIB]) . 
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Note that, as we mentioned in [27*], the parameter sub-space that leads to solutions is neither 
compact nor uniform. The question if it consists from continuous areas or from independent points 
does not have a clear answer [33]. Now, inside these percentages only a small fraction corresponds 
to solutions with < 0, i.e. to oscillatory ones (~ 2% in case IIII Al and ~ 6% in case HUB]) . 
The reason of the significantly smaller appearance of < solutions is that they posses many 
singularities, in comparison with the > case, and therefore it is harder to find a solution with 
no singularity at [0,1]. In conclusion, in total only about ~ 10~^% of the random and uniform 
parameter choices correspond to oscillatory behavior of the metric of the visible brane. In these 
cases, (I2ip and (]26p are numerically verified and the analysis of the cyclic solution is valid. 

Since we have confirmed our analytical calculations numerically, we can proceed to the numerical 
investigation of the general stationary case, i.e. with the full ^ and ([9]) potentials, where analytical 
solutions cannot be obtained. Again, we randomly choose the values of the 8 model parameters 
(10^ octads of Mq, Xq, ao, Ml, Xi, ai, m and A), from a uniform distribution, and we solve the 
transcendental equation system of (|14p and ([6]). Note that the parameter space hyper-cube is taken 
large and symmetric (its edge extends from -10^ to 10^) and we do not impose any constraints on 
the parameter values (assuming for example positive "masses" or opposite tension branes) in order 
to remain as general as possible. In this case only ~ 0.6% of the parameter multiplets corresponds 
to physically accepted solutions, and within them only ~ 4% corresponds to < 0, i.e. in total 
only ~ 10~^% of the parameter choices lead to cyclic branes. We mention that these percentages 
increase radically if we restrict the investigation in specific parameter signs. In fig. [1] we present 
the evolution of the scale factor of the visible brane for one such solution. We conclude that even 
in the absence of analytical calculations in the general-potential stationary case, we do obtain fixed 
and cyclic 3+1 branes and relations (|2ip (or (j26p ) and (|27p . (|28p are satisfied. In other words, if the 
requirements for a complete 5D solution to exist are fulfilled, then (f2T]l (or ([26]) ) and (f27p . (i28|) are 
valid analytically, with the numerics only determining -B(O) and 9. The only numerical restriction 
arises from the presence of B{0) in an exponential, which prevents us from handling arbitrary 
large values as in the analytical calculations. However, this can be solved by the additional scaling 
transformation proposed in [2^. 

Let us make a comment here concerning the "singular" points where the scale factor vanishes. 
These points correspond to the so-called "bounce" , which is always present in all cyclic-cosmology 
models. Although there have been many attempts in the literature to avoid the involved singularity, 



none of them is completed 
quantum fiuctuations 



satisfactory up to now. These approaches, such is the insertion of 



341 ] or the use of loop quantum gravity modifications [35|, could be 
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FIG. 1: Visible brane scale factor evolution for one random solution of the stationary case with the general 
bulk and brane potentials of 0^ and The edge of the parameter space hyper-cube extends from -10'^ to 
10^. The solution arises from m = 81.03, A = -505.8, Aq = -212.5, Ai = -464.2, Mq = 77.94, Mi = 180.7, 
<To = 1.308, (Ti = —3.556 (for simplicity we provide the first four relevant digits only). The obtained B(0) 
and 9 values are correspondingly 9.148 and 0.1522, and t is calculated as e^^^H. All quantities are expressed 
in M5 units. 



included in our analysis, leading to a smoothing out of the behavior of fig. [T] and of relation (j21|) 
(or (j26p ). i.e. making the scale factor non-zero at the bounce. However, in this work we desire to 
present the basic characteristics of cyclic behavior in general collisionless braneworld models, and 
thus we will not examine in detail the (in any case non-complete) handling of the singularity. Our 
model shares this disadvantage of all cyclic models. We will return to this subject in the discussion 
section. 

In the above investigation we have been restricted to stationary solutions of the form (jl2p . case 
in which stable time evolution is implied easily. The question is what can be said about the full 
dynamics of equations ^ and ([5]), where Hq and physical brane distance D{t) can be varying. 



Fortunately, so. 
between them 



utions of the full dynamics seem to consist of stationary ones and the transitions 



2Q|, 



33l |. Therefore our stationary investigation is sufficient. We will return to this 



subject in the next section. 
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V. DISCUSSION-CONCLUSIONS 

In the aforementioned analysis we considered general braneworld models characterized by the 
action ([T]), the conformal metric ([2]), and the general potentials ([8]) and ([9]). Performing both 
analytical and numerical calculations we showed that the full 5D dynamics allows for stationary 
solutions corresponding to oscillatory scale factor of the physical brane and therefore to cyclic 
universes. In statistical terms cyclicity corresponds to ~ 4% of the physical solutions. Our inves- 
tigation is completely 5D, cyclic behavior arises naturally and is induced on the brane by the full 
dynamics, and it is not a result of a modified 4D dynamics, with fine-tuned parameters or specific 
assumptions in the Friedmann equation. Furthermore, we do not use an explicit brane state equa- 
tion, considering just the bulk scalar field (the decays and interactions of which will eventually fill 
the physical brane with the conventional content 136]). As we mentioned in the introduction this 



la]. Indeed, 



full 5D approach is necessary in order to confront the arguments of the authors of 
their allegations that one cannot transit to an effective AD theory (integrating the action over y), 
solve the equations there and then return naively to the 5D description (adding time-dependence 
by hand), are correct. Doing so, the results are not self-consistent (especially the boundary condi- 
tions are not satisfied) and the authors of use this fact as a central argument against the cyclic 
scenario. However, our consistent 51? analysis reveals that cyclic behavior is possible. 

Another important feature of the present study is that cyclic universes do not require brane 
collisions. Thus, we avoid the known problems concerning such a description, which force ekpyrotic 



ne 



aranes do not move at all 



2C 



33]). Furthermore, in our 



model to successively more complicated versions. On the contrary, tl 
and the system is stable (stationary solutions are a stable fixed point 
model, expansion and contraction take place in the 3-1-1 branes, and in all 3-1-1 slices in general, 
while the fifth dimension remains unaffected. The 4 spacial dimensions shrink periodically to an 



ID string and re-expand. This is in a radical contrast with the cyclic mode. 



where the extra dimension is the one that gets contracted (the fifth in [?] or the eleventh in 
Cyclicity seems to re-obtain its "physical" meaning. 

Our 5D investigation is general and does not involve extra assumptions, fine-tunings or specific 
potential forms. We result to periodic, cyclic, homogenous and isotropic universes, where the 
scale factor changes smoothly from expanding to contracting. An observer on the physical brane 
feels successively accelerated expansion, decelerated expansion, turnabout, accelerated contraction, 
decelerated contraction, bounce e.t.c, and a promising signature of the cyclic behavior would be 
the measure of the varying rate of the Hubble constant. The cycles period, given in (j28p . can 



s with extra dimensions, 
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be arbitrary, depending on B(0), i.e. on the value of the warp factor on the physical brane {0 is 
bounded from above and therefore cannot act as a period-decreasing factor). A very interesting 
conclusion comes from the insertion of observational results in our model, which was not made above 
in order to remain as general as possible. Explicitly, if we use the fact that the cosmological constant 
of our Universe is extremely small (~ 0(10"^"^) GeV^), and assuming a reasonable M5 value of 
0(10^^) GeV, the first two boundary conditions in relation ([6]) provide in gene ral a huge value for 



20( 1 ■ which allows us, in 



gB(o) 0(10^^)). This is in consistency with the scaling transformation of 
a solution, to scale the parameters by and add to the warp factor the constant S, and acquire 
another solution. Therefore, the extremely small cosmological constant of the observable universe 
leads the cyclicity period to be around T ~ 0(10^^) years and the maximum scale factor value, 
given by (j27p . to be Umax ~ 0(10^^) m (where the decimal exponents in these rough estimations 
can vary by 1 or 2, depending on -B'(O) and 9 values). Luckily enough, the smallness of the 
cosmological constant excludes oscillatory models with small periods in astronomical terms. In 
more foundational words, the reason that made the cosmological constant that small, is the same 
that makes the cycle period and the size of the Universe that large. 

In this work we have been restricted to stationary solutions, where the subclass of them that 
possesses < corresponds to eternal cyclic behavior with constant period. Numerical investi- 
gation of the full dynamics seem to consist of such stationary solutions and the transitions between 
them 20, [33] . In such transitions Hq on the physical brane can chance sign, leading to a form 



of "chaotic cyclicity", where large intervals of (non-periodic in general) oscillatory behavior could 
be followed by large intervals of conventional evolution and vice versa. In this case, an initial Big 
Bang and/or a final Big Rip or Big Crunch (in conventional terms) could be possible. Another 
interesting possibility would be the exploration of our model with cosmological constants being 
piecewise constant functions of time, reflecting cosmological phase transitions, which could also 
lead to chaotic cyclicity. Note however that numerical confirmation of such behaviors is very hard 
due to the small probability of cyclic stationary solutions (~ 10~^% as we have already mentioned). 
These subjects are under investigation. 

In order for a model to serve as a description of nature, it has to explain the basic physical 
key issues. Especially for cyclic cosmology, amongst others these are the entropy evolution and, 
probably the most pressing issue, that of a fuller understanding of the bounce and the handling of 
the singularity. Our model provides a consistent background for cyclicity and it reveals how such 
a behavior arises from the full 5D dynamics. However, since braneworlds and brane cosmology in 
general arise as limits of a multi-dimensional theory unknown up to now, the 5D results have a 
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phenomenological character and must be considered from this point of view. Definitely, a complete 
explanation and apprehension, and a successful confrontation of the aforementioned subjects, can 
only come through a higher-dimensional, fundamental theory of nature. For the moment we have 
to rely on the relevant research on cyclic cosmology, linearized gravity, M-theory and strings, which 
has improved our knowledge on these issues. These results can be embodied in our analysis. The 
most hopeful effort is the use of quantum fluctuations in order to tame the singularity, which 
effectively is translated into a modification of gravity by the scalar field 



151, 



34( 1 . Alternatively, 



using loop cnantun, g.™ty wecould .odi^- „o„-pertu.bat™ly .he dynan^ical equation, W.g to 

a singularity resolution as in [35] . Concerning the entropy, we could include the relevant discussion 

in our investigation. The argument of the authors of [ll|, [l^ about maximum amount of entropy 

possible in de Sitter spacetime, may lead our model to have a maximum cycle number between 

I I 

10^'' and 10'^'^. However, the idea of the causal patch [12] is probably the best way of handling the 
entropy problem so forth, and there are some interesting recent works on the subject which give a 
boost on cyclic cosmology js?]. 

Let us close this discussion section with some comments on the role of the brane tensions and of 
the bulk cosmological constant in our model. As can be numerically confirmed, setting them to zero 
makes it almost impossible to satisfy the boundary conditions obtaining < and singularity 
absence in [0,1] (this can be achieved only through a careful fine-tuning since our random choice 
procedure gives an one-digit number of such solutions in 10^ parameter multiplets). On the other 
hand, as we showed in IIIIB( in the case where A, Aq and Ai are the only non-zero parameters, an 
~ 10~^% of the random parameter choices, or « 6% of the solutions, correspond to < 0. In 
mathematical terms. A, Aq and Ai are requisite in order to acquire a solution with < in the 
full dynamics, in a natural and not in a fine-tuning way. In terms of physics, it is the dark energy 
that lies in the background of the oscillatory mechanism and allows for cyclicity to realize. Adding 
the fact that it determines the cycles period and the maximum scale factor value, we conclude that 
dark energy is crucial in the described model. This brings it closer to the ekpyrotic paradigm of 
the literature. 

In this work we examine general braneworld models and we show that cyclic behavior can 
naturally arise from the full 5D dynamics. One important feature is that brane collisions are 
not required, on the contrary the branes remain stable, and the cyclicity takes place on the 4Z) 
geometry not on the extra dimension. Another significant result is that the smallness of the 
cosmological constant of the observable universe pushes the cyclic period and the scale factor to 
astronomical large values, an essential requirement for the establishment of cyclic cosmology as a 
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realistic alternative paradigm. Furthermore, we indicate the possibility of a "chaotic cyclicity", 
that is extremely large, non-periodic, cyclic intervals followed by extremely large intervals of 
conventional evolution and vice versa. After these, the model shares both the advantages and 
disadvantages of cyclic cosmology. 
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